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Suppose that a point-like steady source at x = 0 injects particles into a half-infinite line. The 
particles diffuse and die. At long times a non-equilibrium steady state sets in, and we assume that 
it involves many particles. If the particles are non-interacting, their total number N in the steady 
state is Poisson-distributed with mean N predicted from a deterministic reaction-diffusion equation. 

Here we determine the most likely density history of this driven system conditional on observing a 
given N. We also consider two prototypical examples of interacting diffusing particles: (i) a family 
of mortal diffusive lattice gases with constant diffusivity (as illustrated by the simple symmetric 
exclusion process with mortal particles), and (ii) random walkers that can annihilate in pairs. In 
both examples we calculate the variances of the (non-Poissonian) stationary distributions of N. 

PACS numbers: 05.40.-a, 02.50.-r 


I. INTRODUCTION 

Fluctuations of non-equilibrium steady states of driven diffusive lattice gases have attracted a lot of attention in the 
last two decades [J-Q. Although many of these studies assumed purely diffusive, particle conserving dynamics, lattice 
gas models with dissipation have also been investigated wa 3 . The absence of detailed balance makes dissipative 
systems more difficult to handle. There is, however, a steady interest in dissipative models, mostly because of their 
relevance to experiment in such diverse areas as fluid turbulence [l2j , granular gases fl3l f~Mf and many other non¬ 
equilibrium settings in physics, chemistry, biology and engineering. 

The simplest way to characterize a non-equilibrium system is to study its steady state. To maintain a dissipative 
systems in a steady state one must constantly drive it by injecting energy or mass. The total energy or mass content 
of the driven system fluctuates around the mean, and these fluctuations bear the stamp of the non-equilibrium nature 
of the system. In this work we consider three different driven models of diffusive particles where particles can die: 
either individually or by annihilating in pairs. In each of these models the system is driven by injecting particles into 
a half-infinite straight line from a single point-like steady source. The particles diffuse and die so that, at long times, 
a non-equilibrium steady state sets in. If the particles do not interact, their total number N is Poisson-distributed 
with mean N predicted from a deterministic reaction-diffusion equation for this driven system. Our main interest in 
this case will be to find the optimal (that is, most probable) density history of the system conditional on observing a 
given N. We also consider two types of interacting particle models: (i) diffusive lattice gases of mortal particles with 
constant diffusivity but non-trivial fluctuations (as illustrated, for example, by the simple symmetric exclusion process 
jl] with mortality), and random walkers on a lattice that can annihilate in pairs. In both these cases the statistics of 
N is expected to be non-Poissonian, and we calculate the variances describing typical, Gaussian fluctuations of the 
particle number around the mean. We obtain these results by employing (a dissipative extension of) the Macroscopic 
Fluctuation Theory (MFT): a coarse-grained low-noise large-deviation theory that employs, as a large parameter, the 
typical number of particles in the region of interest, see Ref. [l]| for a recent review. The applicability of the MFT in 
the driven systems, considered in this work, demands N 1, and we will work in the parameter regions where this 
condition is satisfied. 

Here is a plan of the remainder of the paper. Section II starts with a brief exposition of the expected, or average 
behavior of a driven system of non-interacting random walkers or Brownian particles that die individually. Then 
subsection IIB presents the MFT formulation of the problem of particle number statistics for a more general family 
of diffusive lattice gases of particles that die individually. For the non-interacting particles, we are able to solve, in 
the same Subsection IIB, the MFT equations exactly. The solution, via the Hopf-Cole transformation, reproduces 
(the large -N asymptotic of) the expected Poisson distribution of N. In addition, it gives the previously unavailable 
optimal density history of the driven system, conditional on observing a given N. Section III deals with two examples 
of interacting particles. The first of them deals with a family of mortal interacting diffusive lattice gases with constant 
diffusivity but non-trivial fluctuations, as illustrated by the simple symmetric exclusion process with mortal particles. 
The second example involves random walkers that only interact via pair-wise annihilation. In these examples a full 
solution of the MFT problem is presently unavailable, and we only calculate the variances of the respective stationary 
distributions of N. Our main results and their possible extensions are briefly discussed in Section IV. 
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II. MORTAL RANDOM WALKERS 


A. Expected density behavior 


Consider a half-infinite one-dimensional lattice with lattice constant a and suppose that a source of particles at the 
origin, x = 0, sets a constant particle number density n o there. The particles perform random walk at x > 0 with 
diffusivity D and die individually with rate /i. When D pa 2 , there is no difference between the discrete random 
walk and continuous diffusion, and the average particle density p(x, t) is governed by the reaction-diffusion equation 


d t p = —pp + Dd^p, 0 < x < oo. 

At long times the average density profile approaches a steady state, independent of the initial condition: 


p{x) = n$e VT5 X . 

Correspondingly, the average steady-state number of particles in this driven system is 


(1) 


( 2 ) 


N= dx p(x) = — no, 

Jo V M 


( 3 ) 


and we assume that this number is much larger than unity. The actual number of particles N in the steady state 
fluctuates around N because of the shot noises of the diffusion and mortality. We are interested in the steady-state 
probability distribution of N. For non-interacting random walkers that die individually this probability distribution 
can be found exactly, by solving the steady-state master equation for the multi-variate probability distribution of 
observing n± particles on site 1, 712 particles on site 2, etc. The solution has the form of the Poisson product measure 
with space-dependent parameters [for D pa 2 , this measure corresponds to the average density profile ©]. This 
leads to a Poisson distribution of N in the steady state. We will proceed, however, as if we were unaware of these 
exact results, and employ instead the MFT: a coarse-grained low-noise theory mostly based on the strong inequality 
N 1. The purpose is two-fold. First, even for the non-interacting particles, the MFT will give the previously 
unknown optimal density history of the driven system, conditional on observed N. The optimal density history is well 
defined only in the limit of TV 1, where it is much more likely than other histories leading to the same N. Second, 
our main interest is in interacting particle models, where exact microscopic results are usually unavailable. In the 
next subsection we briefly discuss the basics of the MFT, and formulate the MFT problem for the particle number 
statistics in a broader context of a family of driven lattice gases of mortal particles. 


B. Macroscopic Fluctuation Theory (MFT) of the Particle Number Statistics 

1. Governing Equations 

When (i) the length scale of interest is much larger than the lattice constant a, (ii) the time scale of interest is 
much larger than the inverse rates of the microscopic processes of diffusion and death, and (iii) the typical number of 
involved particles is much larger than unity, the statistics of large deviations is captured by the MFT jaTTol. Il5l ITU ITiT . 
The starting point of the derivation of the MFT for diffusing and reacting particles is the exact master equation for 
the multi-variate probability distribution of observing a certain number of particles on each cite. Here one can either 
work directly in the physical space, or employ the multi-site probability generating function (that, in the spatially- 
continuous limit, becomes a probability generating functional). Going over to a path-integral formulation, one then 
makes a low-noise approximation by evaluating the path integral by the Laplace method that employs the number of 
particles in the relevant region of space as a large parameter. This procedure yields saddle-point equations (partial 
differential equations) that can be written in a Hamiltonian form: for the density field q{x, t) and a conjugate field 
p{x, t) that plays a role of the “momentum density”. At a qualitative level, the conjugate field p(x, t) describes the 
magnitude of fluctuations. 

If the calculations are performed in the physical space, the saddle point equations (presented here in a form, suitable 
for a class of diffusive lattice gases of particles that die individually), take the form jsl-fTol. ITsl . IToj 

d t q = -pqe~ p + d x [D{q)d x q - a(q)d x p\ , 

dtp = -p{e~ p - 1) - D{q)d x p - ^ cr'(q)(d x p) 2 , 


( 4 ) 

( 5 ) 
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where D(q) is the gas diffusivity, u(q) is (twice) the mobility !lj, and the prime denotes the derivative with respect 
to the argument. The terms proportional to p describe the on-site particle death and its fluctuations. The rest of 
terms describe diffusive transport and its fluctuations. Equations (J4]) and © are indeed Hamiltonian, as they can be 
written in terms of variational derivatives: 


d t q = SH/5p , d t p = —SH/6q, (6) 

where 

POO 

H{q(x,t),p(x,t)} = / dx'H (7) 

Jo 

is the Hamiltonian, and 

= W (e _p - l) - D(q)d x qd x p+^a(q)(d x p ) 2 (8) 

is the Hamiltonian density. Going back to the non-interacting mortal random walkers, we put D(p) = D = const and 
u{p) = 2 Dp, see e.g. Ref. [l]. Then Eqs. © and © become 

d t q = ~pqe~ p + Dd x ( d x q - 2 qd x p) , (9) 

dtp = ~p(e~ p - 1) - Dd 2 p - D(d x p) 2 , (10) 

where 0 < x < oo. The boundary conditions at the particle source are q(x = 0, t) = no and p(x = 0, t) = 0 [Hj. The 
latter condition is quite intuitive: as we demand a fixed (deterministic) value of the density at x = 0, p must vanish 
there. Far away from the source there are no particles. This brings the boundary condition q(x = oo ,t) = 0. 

Being interested in steady state fluctuations, we can assume that, at t = —oo, the system is at its deterministic 
steady state [H, Ho[: q(x,t = —oo) = p(x ), see Eq. ©. We condition the process on observing N particles at some 
finite moment of time that, without loss of generality, we can set to zero. This imposes an integral constraint on the 
solution at t = 0: 

dx q(x,t = 0) = N. (11) 

Analogous integral constraints appear in the MFT formulations of the problem of statistics of integrated current in 
an infinite setting |21] and statistics of particle absorption by an absorber at x = 0 [22j . To account for the integral 
constraint, we should introduce a Lagrange multiplier A and minimize the extended action that incorporates the 
integral constraint. Similarly to Refs. [21] and [22], the action minimization does not change the “bulk” equations 
® and © [or © and mi but yields an additional boundary condition for p(x, t) at t = 0: 

p(x,t = 0) = \9(x), (12) 

where 9(x) is the Heaviside step function, and A is ultimately set by Eq. (fill) [2TI [22t . 

Note that p(x, t) = 0 is an invariant manifold of Eqs. © and ©. The dynamics on this manifold is described by the 
deterministic Eq. ©. This is the relaxation path of the system; it solves the problem in the particular case N = N. 
For N ^ N , the solution of the MFT equations describes the optimal activation path: the most likely density history 
of the driven system conditional on observing N particles. Here p{x,t) ^ 0. Once q(x,t) and p(x,t) are found, we 
can evaluate the action S that yields V(N) up to a pre-exponential factor: 

/ 0 roo 

dt / dx ( pd t q — H) 

-00 J 0 

dx [Dq{d x p) 2 + pq (l - e~ p - pe~ p )] . (13) 




The first term of the integrand comes from the shot noise of diffusion, the second term comes from the shot noise of 
mortality. Rescaling time pt —> t, the coordinate \Jp/D x —> x and the density q/no —> q, one can see that — In V{N) 
obeys a simple scaling relation 


-In P(N) = N f 



(14) 


where f(z) is the large deviation function of the number of particles. Note that no only enters this scaling relation 
through N. 
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2. Particle number statistics and optimal path 


We note that Eq. m is decoupled from Eq. ©■ This decoupling only occurs for non-interacting particles, and it 
gre atly simplifies the problem. Let us perform the Hopf-Cole transformation by introducing Q = qe~ p and P = e p — 1 
[18j . The generating functional of this canonical transformation can be chosen to be 


F{q(x,t),Q{x,t)} 





(15) 


In the new variables Q and P the Hamiltonian is H{Q(x,t), P(x,t)} = J dx'H , where 


U = -HQP- Dd x Qd x P. 

As a result, the MFT equations become linear and fully decoupled: 

d t Q = -pQ + DdlQ, 
d t P = pP-Dd 2 x P. 


(16) 

(17) 


Note that these equations for Q and P arise immediately, when one employs the Laplace method for the evaluation 
of the path integral in the formalism of probability generating functional |18| . 

In the new variables Q and P, the boundary and initial conditions are: 



Q(0, t) = no and Q{x, — oo) = p(x) 

(18) 

for Q, and 

P(0, t) = 0 and P(x, 0) = (e A — l) 6(x) 

(19) 

for P. As a 

result, Q(x, t) is invariant in time, 



Q(x, —oo < t < 0) = p(x), 

(20) 

while 

P(x, oo < t < 0) — (e x l) erf ^ . 

(21) 


Now we can determine the optimal path in the original variable q: 


q(x, t ) = Q{x)[ 1 + P(x, t)] = n 0 e~ 
Using Eq. (fill) , we find A = In (N/N), so 


1 + (e x - 1) e pt erf 




( 22 ) 


q{x,t) =n 0 e V S 


l + (^-l 


N 


z pt erf 


yf^PDi) _ ' 


(23) 


It is easier to calculate the action in the new variables Q and P where, as one can show by a direct calculation 
the action is equal to the increment of generating functional F from Eq. (USD: 


S = F{q(x , 0), Q(x, 0)} - F{q(x, —oo), Q(x, —oo)}. 


(24) 


After some algebra, this gives 

-In V{N) ~ S 



q(x ,0)In 


g(jE.O) 

Q(x, o) 


N 

N ln^ 
N 


-N + N, 


q{x,0) + Q(x, 0) 


(25) 


where N is given by Eq. ©. That is, the large deviation function f(z) from Eq. (fLU) is equal to f(z ) = z In z — z + 1. 
The distribution (l?5l) coincides with the N 1, N 1 asymptotic of the Poisson distribution with mean N, as to 
be expected. In particular, the variance of this distribution coincides with the mean: 


U RW =N = 



(26) 
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FIG. 1: (Color online) The optimal density history for N/N = 3 (the left panel) and N/N = 1/3 (the right panel). The times, 
rescaled by the decay rate p , for both panels, are —oo (solid line), —0.4 (dashed line), —0.1 (dotted line) and 0 (dash-dotted 
line). At t = 0 a density jump develops at the origin, so that the required number of particles N effectively comes from a 
deterministic density profile with a different density at the source. 


Now let us return to the optimal path (E^l) that has been previously unknown. Although the statistics of N is 
time-independent, the optimal path does depend on time. Furthermore, the activation path does not coincide with 
the time-reversed relaxation path, obtained by solving the deterministic reaction-diffusion equation m back in time. 
This is a clear signature of non-equilibrium. Notice also that, in order to ensure an unusually large or small number of 
particles at t = 0, the fluctuations create a boundary layer in the density profile at the particle source. This boundary 
layer becomes a density jump at t = 0, 


q(x>0,t = 0) = r^p(x), (27) 

so that the effective boundary condition is q(x —> Off = 0) = uqN/N , whereas the bulk of the gas particles behaves 
deterministically. These features can be seen on Figure Q] which shows the optimal density histories described by 
Eq. ([25)1 . The left and right panels correspond to N = 3 N and N = N/ 3, respectively. These results are both 
unexpected and instructive. 


III. INTERACTING LATTICE GASES 

Now let us consider interacting lattice gases of mortal particles. They may have different a(q ), but for simplicity we 
will continue to assume a constant diffusivity D = const. For such gases the deterministic equations ©-© continue 
to hold, while the MFT equations read 

d t q = —pqe~ p + d x [Dd x q - a(q)d x p] , ( 28 ) 

d t p = -p(e~ p - 1 ) - Dd 2 x p- ^a'(q)(d x p) 2 , ( 29 ) 

with the same boundary conditions as before. Once q(x,t) and p(x,t) are known, the probability distribution V(N) 
can be evaluated from 

/ 0 poo 

dt / dx 

-oo J 0 

Rescaling time pt —> t and the coordinate \Jp/D x —> x : we obtain a scaling relation 

-\nV(N) = Nf(^L,n o y ( 31 ) 

where no enters both through N and separately [ 24 | ■ 

It does not seem possible to solve Eqs. (E51) and m and determine the large deviation function f(z , no) analytically 
for a general cr(q). Typical, Gaussian fluctuations of the number of particles around N are given by a quadratic 


-cr(q){d x p) 2 + pq(l-e p - pe p ) 


(30) 
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asymptotic of f(z) at z close to 1. This asymptotic can be found relatively easily via a perturbation theory around 
the deterministic steady-state solution ©. This theory employs the Lagrange multiplier A as a small parameter [23| . 
We set 


q = p(x) + Aqi + A 2 q 2 + ... (32a) 

p = Xpi + A 2 j> 2 T • ■ ■ (32b) 

and plug these expansions into Eqs. K51) and (1251) . The first-order equations are 

(d t + p — Ddl) qi = pp{x)p 1 - d x [<j{p)d x p 1 ], (33a) 

(d t - p + L><9 2 ) Pi = 0- (33b) 

The equation for pi is independent of <r(q), and it is decoupled from the equation for q\. Therefore, we can solve it 
immediately, with the boundary conditions pi(0, t) = 0 and pi(x, 0) = 9(x). The solution is 

Pi(x,t < 0) = e Mt erf . (34) 

Now we could plug this expression in Eq. (133al) and solve for qi with the boundary conditions qi(0,t) = 0 and 
qi(x, — oo) = 0. This is unnecessary, however, for the purpose of computing the variance of V(N), because the latter 
is independent of q\. Indeed, we have 

A 2 /'° r°° 

- In V(N)~S = — / dt dx [a{p){d x p{) 2 + pppl] + 0(A 3 ). (35) 

* ,7-00 Jo 

Correspondingly, the variance is equal to 

/ 0 noo 

dt / dx [cr(p)(d x pi) 2 + ppp'l] , (36) 

-oo J 0 

with p(x) from Eq. (J2J and pi(x,t) from Eq. (l34l) . Equation (l36l) does not demand a knowledge of qi(x,t) and holds 
for any mortal diffusive lattice gas with D = const. 

One well-known example of a gas with constant diffusivity D but non-trivial fluctuations is provided by the SSEP, 
where each particle can randomly hope to a neighboring lattice site if that site is vacant. If it is occupied by another 
particle, the move is forbidden. For the SSEP one has <r(p) = 2Dp(l — p) [l]. Here we set the lattice constant a = 1, 
so that the particle density at the source 0 < n 0 < 1 is dimensionless. Evaluating the double integral in Eq. (1M1) in 
this case (see Appendix [All . we obtain 


Vssep 



(37) 


As Vssep 7^ N, V(N) is non-Poissonian. As expected on the physical grounds, Vssep is smaller than the variance of 
the total number of non-interacting random walkers with the same no, see Eq. (1261) . so the distribution is narrower 
than the Poisson distribution with the same mean. The two variances coincide in the limit of no —> 0, where exclusion 
effects in the SSEP are negligible. That V(N) is non-Poissonian is not surprising, but even its variance has been 
previously unknown. 


IV. ANNIHILATING RANDOM WALKERS 

An annihilating random walker (ARW) is immortal when it is alone, but two ARWs on the same lattice site can 
annihilate, 2 A —> 0. Let a be the annihilation rate constant. When the diffusion is sufficiently fast [see the criterion 
m below] the average particle density p(x, t ) is governed by the continuous reaction-diffusion equation [l9j 

d t p = — ap 2 + Dd x p , 0 < x < oo. (38) 

The particle source at x = 0 fixes a particle density n 0 of the ARWs there. With this boundary condition, the 
steady-state average density profile is 


ano 

6 15 


p{x) = n 0 1 + 


-2 


(39) 
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it falls off much slower than the exponential profile fl2J. For the continuous reaction-diffusion equation to be valid, it 
is necessary that the characteristic length scale ~ (D /aji-o) 1 ^ be much larger than the lattice constant a: 


D 

ano 


» a. 


The average steady-state number of particles is 


1/2 


(40) 


(41) 


and we assume N 1. The MFT equations for this system can be safely derived from the exact master equation 
for the multi-variate probability distribution by assuming that the typical number of particles on each lattice site is 
much larger than unity, leading to the strong inequality no a ~S> 1 [SGI. We believe, however, that it is actually 
sufficient to require a weaker condition N 1, alongside with the condition (HUl) . The MFT equations are IS 


whereas 


is the Hamiltonian density G1G1- The first term comes from the on-site annihilations, the second and third terms 
come from the diffusion. Correspondingly, 


d t q = -aq 2 e 2p + Dd x ( d x q - 2qd x p ), 

(42) 

d t p = - aq(e~ 2p - 1) - Dd^p - D(d x pf, 

(43) 

'H(q,p) = T^uq 2 (e _2p - l) - Dd x qd x p + Dq(d x p) 2 

(44) 


- In■p(iV) ~ S' = 


dt 


dx 


Dq(d x p) 2 + £ q 2 (1 - e~ 2p - 2pe~ 2p ) 


As one can check, by performing rescalings described by Eq. (l48l) below and additional rescaling q/no —> q, 

- In P(N) = Nf 


S) 


(45) 


(46) 


Here no only enters through N, as for the random walkers who die individually, cf. Eq. m 

As in Sec. nm we can calculate analytically the variance of the total number of ARWs in the steady state. We 
make the ansatz (13 2 all and (132bl) in Eqs. (l42l) and (143II and obtain, in the first order in A <C 1, 

(47a) 
(47b) 


[d t + 2 ap(x) - Dd%] q x = 2ap 2 (x)p 1 - 2Dd x {p{x)d x pi\, 
[d t - 2ap{x) + Ddl] pi = 0. 


As for the mortal SSEP, Eq. (I47bl) for pi is decoupled from that for q lt and its solution suffices for computing the 
variance we are after. Let us reverse and rescale time and transform the coordinate: 


anot 


i i , an ° 

and y = 1 + \ —— x, 


6 D 


so that Eq. (I47bl) becomes 


a , 12 P 1 02 

d T p i 4 - 5 - = d p 1 . 

y 


(48) 


(49) 


We need to solve it for f < y < 00 and 0 < r < 00 subject to the boundary condition p\{y = l,r) = 0 and initial 
condition 


Pi(y,r = 0) = 9(y- 1). 


(50) 


The rescaled problem for pi{y,r) is parameter-free, and we solve it in Appendix [Bj Once pi(y,r) is found, we can 
calculate the variance of N: 


/ 0 poo 

dt / dx [2Dp(d x pi) 2 + 2ap 2 pl\ 
-00 J 0 


Jo J 1 y L 


(a \2 . 

(d y p 1 + — 

y 2 


( 51 ) 
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As expected from Eq. (l46l) . Varw is proportional to N. Using Eq. (IB 101) of Appendix iBl for p(y, r), we can represent 
the proportionality coefficient (a dimensionless number of order unity) as a quadruple integral. The integration over 
r is elementary. The integration over y is very tedious, but can be performed explicitly with “Mathematica”. We 
evaluated the remaining double integral numerically, leading to Varw — 0.78 N. 

As a check, we also solved Eq. (l49l) numerically in the region 1 < y < L and 0 < r < T with the boundary 
conditions p\{y = l,t) = 0 and d v p\{y = L,t) = 0 and initial condition ijhUjl . taking L and T sufficiently large. Then 
we used the numerical solution to compute the double integral in the second line of Eq. m- The result comes quite 
close, Varw — 0.77 N. As Varw < A) the distribution V(N) is narrower than a Poisson distribution with the same 
mean. 


V. SUMMARY AND DISCUSSION 

This work addressed the statistics of the total number of particles N that are present at any chosen time in the 
steady state of a driven lattice gas composed of mortal diffusing particles. The formalism we used is that of the 
dissipative Macroscopic Fluctuation Theory (MFT). For non-interacting random walkers who die individually, the 
MFT formulation of the problem is exactly soluble and yields the expected Poissonian statistics of N with the mean 
predicted by the simple reaction-diffusion equation ©• It also provides a fascinating and instructive visualization of 
large deviations of N in the form of the optimal density history of the driven system conditional on N. For interacting 
diffusing particles we calculated the variance of the distribution of N, and found that the distribution is narrower 
than a Poisson distribution with the same mean. 

The variance calculations that we showed here is a first step towards studying the complete statistics of N. Extending 
our perturbation theory for the MFT to higher orders in A, one should be able to compute several higher moments 
of V{N), as it has been recently done in the problem of melting of an Ising quadrant [25}. We also note that it is 
possible to compute the distribution of N numerically by solving the full MFT equations with the proper boundary 
with the Chernykh-Stepanov iteration algorithm |26|. This algorithm was originally developed for evaluating the 
probability distribution of large negative velocity gradients in the Burgers turbulence. Later on it was used in studies 
of different types of large deviations in diffusive lattice gases, with and without on-site reactions fl8l. Hfil. |23[ l27l - [30} . 
This algorithm is much more computationally efficient than microscopic stochastic simulations. 

It would be very interesting, and challenging, to directly probe the tails of V(N), that are beyond the reach of the 
small-A perturbation theory. For the SSEP involving immortal particles, y, = 0, the limit of very large transferred 
mass, in an infinite system, can be described by neglecting the term —Dd x qd x p in the Hamiltonian density ©. The 
ensuing reduced MFT equations turn out to be exactly soluble [2)1 [31]. Whether a similar reduction is possible in 
the problem of extreme statistics of the total number of interacting mortal particles is an open question. 

On a more general note, understanding non-equilibrium systems requires, among other things, intuition which one 
acquires by learning from examples. The prototypical dissipative systems, considered in this work, are helpful in 
gaining such an intuition. 
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Appendix A: Calculating the variance for the SSEP 

Here we evaluate the integral in Eq. (IMl) for the SSEP. We start with calculating the integral 

pO poo 

h= dt dx a(p)(d x pi) 2 

' —oo JO 


J dxj dt2Dnoe ^1 — no e ^ 

-^ x (l - no e-^ x ) K 0 (, 


~7Bt + 2 f J ' t 


7T Dt 


4 n 0 


7r 


(Al) 
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where Kq(- ■ •) is the modified Bessel function of the second kind. Evaluating the remaining integral, we obtain 


r D ( 4-\/3 2 n 0 


Now we evaluate the integral 


h = 


pOO 

dt / dx npp\ 

Jo 


— oo JO 

r° r 
dt / 

— OO J 0 


dxiinoe-V^ x+2 ^ erf 2 



4 D 

n 0 
A* Jo 


V^Dt 

/*oo /»oo _ 

/ dz\fze~ 2z / du e~ 2 ^ zu erf 2 it. 

Jo Jo 


With some patience, this double integral can be evaluated using the book of integrals [32|, and the result is 


(A2) 


(A3) 


(A4) 


Summing up I\ and I 2 we obtain Eq. (ED. 


Appendix B: Finding pi(y,r) for annihilating random walkers 


To find the spectrum of the problem and the proper eigenfunctions, we make the ansatz Pi(y,r) = tl>(y,T)e r 
and arrive at the equation 


d 2 ip(y,T) 


dy 


12 


2 +( r -^)V>(?/,r) = o 


(Bl) 


that we need to solve with the boundary condition tjj{y = 1,T) = 0. Equation (IBID is the Shrodinger equation for a 
quantum particle with energy T 2 in the potential 


v(y) = \%’ y>l : 

00 , y < 1. 


(B2) 


The spectrum of the problem is continuous, 0 < T < 00 . Two linearly independent solutions of Eq. m can be 
chosen as 


and 


i’Jy, r) = 


V >2 (y,r) = 


Ty (T 2 y 2 — 15) sin(T?/) + 3 (2 T 2 y 2 — 5) cos(Ty) 
T 3 y 3 

Ty (r 2 y 2 — 15) cos(Ty) — 3 (2 T 2 y 2 — 5) sin(Ty) 

fv ' 


(B3) 


(B4) 


At T —> 00 (or at y —> 00 ) these solutions become sin(Ty) and cos(Ty) as expected. All the eigenfunctions, vanishing 
at y = 1, can be written as 


where 


4>{y, r) = a(T)<j}(y, r), 


4>{yX) = ^i(2/,r)^ 2 (i,r) - i/>i(i,r)^ 2 (i/,r) 


(B5) 


(B6) 


and a(T) is a yet undetermined amplitude. The eigenfunctions (IB5D are orthogonal, and they can be normalized as 
follows: 


/ OO 

dy^(y,T)^(y,T')=S(T-T'), 


(B7) 
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where S is Dirac’s delta function. As a result, 

pOO pOO 

a~ 2 ( T)= dV dy<t>(v,T)<t>(vX)- 

Jo j 1 

Evaluating this double integral with a help of “Mathematica”, we obtain 


(B8) 


a(r) 


The solution for pi(y,r) can be written as 


[2 r 3 

V 71- Vr 6 + 6r 4 + 45T 2 + 225' 


(B9) 


poo 

Pi(y, T )= dTA(T)'i/j(y,T)e~ r2T , (BIO) 

Jo 

where A(T) is the projection of the initial condition (f50l) on the normalized eigenfunctions (|B5p . That is, A(T) = 
dy^(y,T). This integral can be also evaluated with “Mathematica”, resulting in a tedious formula 

4F 3 4- 90r + 6Ci(r) [3 (5 — 2r 2 ) sinT + r (r 2 - is) cosF] — 3r (r 2 — 15) [tt — 2Si(r)] sinT — 9 (2r 2 — 5) [n — 2Si(r)] cosT 

2v / 27rr Vr e + 6r 4 + 45r 2 + 225 


where 

Si(*) = f'^dy and Ci(*) = - f = dy 
Jo y Jz y 

are the sine and cosine integrals, respectively. Now pi(y,r) in Eq. (IB10I) is fully determined in terms of a single 
integral over T. 
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